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Abstract: We study the classical version of the 120° -model. This is an attractive nearest-neighbor 
system in three dimensions with XY (rotor) spins and interaction such that only a particular pro- 
jection of the spins gets coupled in each coordinate direction. Although the Hamiltonian has only 
discrete symmetries, it turns out that every constant field is a ground state. Employing a combi- 
nation of spin-wave and contour arguments we establish the existence of long-range order at low 
temperatures. This suggests a mechanism for a type of ordering in certain models of transition- 
metal compounds where the very existence of long-range order has heretofore been a matter of 
some controversy. 
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1 . Introduction 

1.1 Overview. 

For attractive classical spin systems with ground states related by an internal symmetry, ordering 
usually occurs by one of two mechanisms: The existence of surface tension between thermal 
perturbations of the ground states or condensation of spin-wave deviations away from the ground 
states. The former is most common in models where the internal symmetry is discrete, while 
the latter circumstances are best exhibited in systems with continuous symmetries. This paper 
will be concerned with an attractive spin system — the so called 120°-model — which displays 
characteristics reminiscent of both phenotypes. A related model of this sort — the so called orbital 
compass model — will be the subject of a continuation of this paper [4]. A common feature of both 
systems is that the presence/absence of long-range order is all but readily apparent. 

To underscore the above (admittedly vague) allegations, let us introduce the formal Hamilton- 
ian of the 120° -model: 

^ = - sfUf + - + [SP - Sl%f}. (1.1) 

r 

Here r is a site on the cubic lattice Z^, the S/s wee the usual XY-spins, namely two-dimensional 
vectors of unit length, and e^;, e^ and are the lattice unit vectors in the three coordinate direc- 
tions. To define the quantities si-^^ , Sr'"^ and Sr'^^ , let a, b and c denote three vectors on the unit 

circle evenly spaced by 120°. Then Sr^^ = Sr a. and similarly for Sr^^ and Sr'^\ We have J > 
so the interaction is ferromagnetic. 

As is manifestly obvious from (1.1), any constant spin field is a ground state and since we are 
dealing with continuous spins, no contour-based argument readily suggests itself. (As we shall 
see later, there are also other ground states, but these need not concern us at the moment.) On 
the other hand, due to the directional bias of the coupling, a naive spin-wave argument based 
on the use of infrared bounds [20, 22, 23, 25] results in divergent momentum-space integrals. In 
particular, as we later show, the spherical version of this model has a free energy that is analytic 
at all temperatures. Worse yet, the rigorous version of a disorder-by-spin-wave argument, the 
Mermin-Wagner theorem, requires the continuous symmetry to be present at the level of the 
Hamiltonian, which here is simply not the case. Thus, the system in (1.1) is right on the margin. 

The main goal of this paper will be to establish long-range order in this model. (Precise defi- 
nitions will appear at the end of this section; precise statements of the theorems will appear in the 
next section.) The mechanism for ordering involves the combination of different aspects taken 
from both of the classic types of arguments. Specifically, on the basis of a realistic spin-wave 
calculation we show that, for all intents and purposes, most of the ground states are destabihzed, 
leaving us with only a manageable number of contenders. Among the survivors, a surface ten- 
sion (with some unusual features) is established. Thereafter, via arguments which are relatively 
standard, the existence of multiple states at low temperatures can be concluded. 

The described reduction of the ground state degeneracy by accounting for the free energy of 
the excitations is reminiscent of the problems analyzed previously in [7, 16, 29]. (In our cases. 
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the role of excitations is taken by the spin waves.) However, the model studied in this paper 
presents us with several novel features. For instance, unlike in [7, 16, 29] which focused on 
discrete spin systems with "stratified" ground states, here we are dealing with a continuum of 
homogeneous ground states related by a continuous internal symmetry. Incidentally, "stratified" 
ground states also exist in our systems, see Sect. 1.3. Here these must be ruled out on the basis of 
a modified spin-wave calculation which accounts for the free energy carried by deviations from 
inhomogeneous background. 

Although the authors would have been proud to stake the claim of having concocted a model 
system with such an esoteric mechanism of ordering, it turns out that interest in the 120° -model — 
as well as the closely related orbital compass model — is not entirely academic. Indeed, both 
systems arise naturally in the study of transition-metal compounds. Here magnetic order of some 
type has been firmly established by experimental methods, but the nature and the mechanism 
for the order is unclear. The problem persists up to the theoretical level; the question whether 
any interacting model based on the physics of transition-metal orbitals is capable of supporting 
long-range order has heretofore been a matter of controversy. From the present paper we now 
know that, at the level of finite-temperature classical spin systems, ordering indeed occurs for 
the 120° -model. This strongly suggests (but of course does not prove) that a similar ordering is 
exhibited in the quantum and itinerant-electron versions of these models. 

The rest of this paper is organized as follows. In Sect. 1.2 we describe the physical origins 
of these problems. A precise definition of the classical 120°-model is given and the ground 
states are discussed in Sect. 1.3. In Sect. 2.1 we state our main result concerning the existence 
of phase transition in the 120° -model while in Sect. 2.2 we outline the principal ideas of the 
proof of long-range order. The actual proofs are given in Sect. 3. The techniques we employ are 
contour methods based on chessboard estimates but the infinite degeneracy of the ground states 
also requires us to perform some intricate spin-wave calculations. These technical details are the 
subject of Sects. 4—6. Sect. 7 collects some observations concerning the spherical version of the 
model at hand. 

1.2 Quantum origins. 

In the standard description of electrons in solids, it is often the case that the accumulation of 
itinerant charges is heavily disfavored. This (presumably) results in localized electrons which 
interact only via spin exchange. In the circumstances which are most often studied, only a sin- 
gle orbital is available at each site, which produces an effective antiferromagnetic interaction. 
However, in transition-metal compounds (e.g., vanadates, manganites, titanates, cuprates, etc.) 
there are multiple essentially-degenerate orbitals any of which could be occupied. In particular, 
if the transition metal ion interacts with a local environment which is of octahedral symmetry, 
the 3d-quintet of the transition-metal ion is split into a low-lying triplet — ^the t2g orbitals — and a 
pair — the Cg orbitals — of considerably higher energy. 

In the absence of any other significant effects, one circumstance which is amenable to further 
approximation is when there is but a single electronic degree of freedom per site. The two ob- 
vious distinguished cases are the eg and t2g compounds. The former will come about under two 
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conditions: First, if the t2g orbitals are filled and there is one extra electron per site to occupy 
the Gg orbitals. Second, same as above but here there are three electrons (out of possible total of 
four) in the eg orbitals and the role of the single electronic degree of freedom is played by the 
hole. The latter cases, the t2g compounds, occur if the eg orbitals are empty and there is a single 
electron or a hole in the t2g orbitals. It appears that the situations leading to e^-type compounds 
are far more prevalent. 

In any of the above circumstances, one can write down the inevitable itinerant electron model 
describing the spins and orbitals. After the standard superexchange calculation/approximation — 
analogous to that which relates the single-orbital models to the Heisenberg antiferromagnets — we 
arrive at a problem which involves "only" quantum spins. Of course, in these models there will be 
two types of quantum states. Namely, those corresponding to the actual (electronic) spin degrees 
of freedom and those corresponding to the occupation numbers of the dynamical orbitals. The 
resulting system is described by the Kugel-Khomskii [34] Hamiltonian 

a r 

Here the interaction takes place at the neighboring sites of the cubic lattice representing the 
positions of the transition-metal ions, the object cr^ is the triple of the usual Pauli matrices acting 
on the spin degrees of freedom at the site r and the tt^ are pseudospin operators acting on the 
orbital degrees of freedom at the site r. As usual, the vectors e^;, Cy and e^ are the unit vectors in 
the principal lattice directions. 

Depending on which of the orbitals play the seminal role, the two choices for the orbital pseu- 
dospins are 

-a V3c7^), if a = x,y, 

< = < , (1.3) 
[^cr^, ifa = z, 

for the Cg-compounds, while 

vr" = ^it", for a = x,y,z, (1.4) 

for the t2g-compounds. The former choice gives rise to the 120° -model and the latter to the orbital 
compass model. 

The question of obvious importance is to prove/disprove the case for ordering of the spins 
or orbitals in these models. In this vain, it should be remarked that the orbital-compass version 
of the Kugel-Khomskii Hamiltonian — if reformulated back in the language of itinerant-electron 
model — ^has some unapparent symmetries. For instance, as pointed out by Harris et at [28], the 
total spin of electrons in a-orbitals at sites of any plane orthogonal to the direction represented 
by a is a conserved quantity. On the basis of these symmetries, a Mermin- Wagner argument has 
been used [28] to show that, in the three-dimensional system, the spin variables represented by Cr 
in (1.2) cannot order. 

Notwithstanding the appeal of this "no-go" result, we note that the absence of spin order does 
not preclude the more interesting possibility of orbital ordering in these systems. Indeed, on 
the experimental/theoretical front, it appears that there is a reasonable consensus "for" orbital 
ordering; the references [9, 11, 12,21,30,39,41] constitute works which support this picture while 
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the references [1, 31, 35,46] offer arguments that dispute or down-play the role of orbital order 
in the magnetic properties of transition-metal compounds. We refer to (slightly biased) review 
articles [8,47] for more information. 

In order to study the phenomenon of orbital ordering in the context of Kugel-Khomskii mod- 
els, the interactions are often further reduced. Neglecting all sorts of terms including all terms 
pertaining to intrinsic spin, the resulting orbital-only model has the Hamiltonian 

■^ = '^EEW+e„- (1-5) 
a r 

Here, as before, tt" are as in (1.3) for the 120° -model and (1.4) for the orbital compass model. 
Full physical justification of these approximations goes beyond the scope of this paper. 

Interestingly enough, the Hamiltonian (1.5) for the 120° -case can be arrived at by entirely 
different means. In particular, among the other "competing" mechanisms so far omitted from the 
discussion is the Jahn-Teller effect which refers to further distortion of octahedral geometry of 
the "crystal field" surrounding the transition-metal ions. On the basis of symmetry considerations 
it has been argued [32] that, in the e^-compounds, this will lead to an effective interaction among 
the nearby orbitals which turns out to be exactly of the type (1.5). In the rare cases of the t2g- 
compounds with Jahn-Teller effects, it turns out that yet another Hamiltonian emerges. In the 
t2g-cases the interplay of the two interactions must be properly accounted for; in contrast to the e^- 
situations where, no matter what, we get the 120° -model. For these and other reasons — the latter 
mostly concerning the "degree" of difficulty — the remainder of this paper will be focused on 
the 120°-model. 

1.3 The classical models. 

The classical versions of the above orbital models can be obtained from their quantum counter- 
parts by replacing the operators tt^ by appropriate projections of the classical spin variables Sr, 
which live on the unit sphere in M". A standard justification for the classical approximation is 
via the "S oo" limit; cf [18,36,45] and also [19,37,38] for some results in this direction. As 
was the case for the quantum systems, there are two major types of models under consideration: 
Classical 120°-model and classical orbital compass model. We proceed with formal definitions. 

Let denote the three-dimensional cubic lattice and let Sr, where r G Z^, be unit vectors 
in M?. We let a, b and c denote three evenly-spaced vectors on the unit circle, for instance, 

a=(l,0), b=(-i,^) and c=(-^,-^), (1.6) 

fa) 

and define the projections Sr = Sr-a, where the dot denotes the usual dot product, and similarly 
for Sr^^ and Sr'^\ In this notation, the (formal) Hamiltonian of the 120° -model is given by 

= -J^^isPs^,^ + Si'^sl%^ + sPsj^J, (1.7) 

r 

with again J > 0. For convenience we will sometimes label the lattice direction and the spin 
direction with the same index; i.e., Sr'^\ a = 1,2, 3, meaning, e.g., Sr^^ for a = 2, etc. Then 
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(1.7) can be written 

j^ = -jysj^''^slll . (1.8) 

We remark in passing that (1.8) is also the form of the orbital compass Hamiltonian, only in this 
case the Sr, r G Z^, are genuine (three-component) Heisenberg spins and the upper index of the 
spin stands for its Cartesian component. 

Remark 1. The 120° -model (as well as its orbital compass counterpart) can be generalized to 
hypercubic lattices in other dimensions as well as to other graphs. For instance, in four spatial 
dimensions the spins are from the unit sphere in and the interaction in the various lattice 
directions is via the projections of the spins onto the vectors pointing from the origin to the 
vertices of an appropriately centered tetrahedron. However, these variant situations are fairly 
difficult geometrically and since they do not always correspond to the structure of the original 
quantum-spin model, we will not consider them in this paper. 

The salient feature of the 120° -model (as well as the orbital compass model) is that the ground- 
state space of the Hamiltonian is infinitely degenerate. This is manifest if we write the Hamilton- 
ian in the form (1.1) which follows immediately from (1.7) by the fact that for any Sr from the 
unit circle in M?, 

[sPf + [sl^^f + [sPf = l. (1.9) 

It is now apparent that any constant vector field receives the minimum possible energy — ^namely 
zero — from the Hamiltonian in (1.1). 

Unfortunately, as we remarked before, the ground state situation is further complicated by the 
fact that the constant configurations are certainly not the only minimum-energy states available 
in this system. For instance, it is easy to verify that, starting from a constant configuration, the 
reflection of all spins in an xy-plane "through" vector c preserves the overall energy. (Here, the 
c-projection is not affected by this procedure and the a and b-projections just swap their roles.) 
Hence, plenty of other ground states can be obtained from the constant ones by reflecting all spins 
in a collection of parallel lattice planes; see Fig. 1 for some examples. Notwithstanding, as will 
be proved later, these non-translation invariant ground states are disfavored by the onslaught of 
positive temperatures. 

Remark 2. The ground state situation is yet more intricate in the orbital compass model which is 
foremost among the reasons that our analysis of this system was postponed. 

1.4 Gibbs measures. 

The (still formal) Hamiltonian in (1.1) can be used to define the Gibbs measures for the corre- 
sponding spin system. Explicitly, let A C be a finite set and let dA denote the set of sites 
in \ A that have an edge with one endpoint in A. Given a spin configuration Sa in A and a 
boundary condition Sqa on dA, we let J^(5A|5aA) be the restriction of the sum in (1.1) to r 
and a such that r G A or r -|- Sq G A (or both). Then the finite- volume Gibbs measure in A with 
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Figure 1. A picture demonstrating the discrete symmetries of the ground states in the 120°- 
model on a cube with one spin fixed. Here the horizontal and vertical directions correspond to 
the y and z-axes, respectively; the front face of the cubes is perpendicular to the s-axis. The 
upper-left cube is simply the homogeneous ground state, the upper-right cube has a spin reflection 
in the a-direction as one moves in the Sa; -direction across the yz-midplane. The bottom cubes 
have analogous b and c-reflections. The structure of any (global) ground state is demonstrated by 
checking for consistency between all neighboring cubes. 



boundary condition S^a is a measure on the configurations Sa = (S'r)reA given by 



n ^i^Sr). 

reA 



(1.10) 



Here Q is the Lebesgue measure on the unit circle (in other words, O, is the a priori spin distribu- 
tion) and Z^^^ is the corresponding normalization constant. 

Regarding these measures as the so-called specifications, the DLR-formalism can be used to 
define the infinite- volume Gibbs measures (aka Gibbs states). ExpUcitly, the latter are probability 
measures on configurations Sr for r G Z'^, whose conditional probability in a finite volume given 
a boundary condition Sq\ is the measure (1.10), for almost every Sq\. We refer to [26] for a 
comprehensive treatment of these concepts. To adhere with mathematical-physics terminology, 
we will denote expectation with respect to the infinite- volume Gibbs measures by (— ). 
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2. Main results 

Here we state the main theorem of this paper and provide heuristic reasoning for the existence of 
long-range order in the system at hand. We also provide some more discussion and remarks on 
literature concerning the related problems that have previously been studied. 

2.1 Long-range order in the 120-degree model. 

Let w^, r = 1, 2, . . . , 6, denote the six vectors on the unit circle in corresponding to the six 
sixth-roots of unity. Explicitly, we define 

w^ = (cos(fr),sin(fr)), r = l,2,...,6. (2.1) 

The principal result of this paper is then as follows: 

Theorem 2.1 Consider the 120° -model with a fixed coupling constant J > 0. Then there exits 
a number Pq G (0, oo) and a fiinction f3 i— > e(/3) G [0, 1) with e(/3) as f3 ^ oo such that 
the following is true: For all (3 > Po there exist six distinct, infinite-volume, translation-invariant 
Gibbs states (— )^ j, with r = 1, 2, . . . , 6, such that 

{^r-Sr)lj>l-e{P), r=l,2,...,6, (2.2) 

is valid for all r € Z"^. 

We note that once (w^ • SrYp ; / 0, we must have that {SrYp j ^- Consequently, (2.2) 
implies the existence of a long-range order because at /3 ^ 1, the standard high-temperature 
expansions (or Dobrushin uniqueness techniques, see [44, Theorem V. 1.3]), imply that = 
in any Gibbs state (— )/3,j. Moreover, as /3 — oo, the measure corresponding to (— )^ j gets 
increasingly concentrated around w,-. 

Theorem 2.1 is proved in Sect. 3.2 subject to some technical claims whose proof is postponed 
to Sect. 6.3. 

2.2 Spin-wave heuristics. 

Here we provide a heuristic outline of the spin-wave reasoning which ultimately leads to the proof 
of the above theorem. The precise version of the argument is given in Sects. 4 and 5. 

The starting point of our analysis differs in perspective from the usual sort of spin-wave ar- 
guments which have previously been the subject of mathematical theorems. In the standard ap- 
proaches, one attempts to rewrite the full Hamiltonian as a "spin-wave" Hamiltonian, carry out 
a calculation and control the errors later (if at all). An extreme example of this is the spherical 
model whose working definition is "the spin system for which the spin-wave approximation is 
exact." However, as alluded to previously, this sort of spin-wave approximation is inadequate to 
capture the essential features of the problem at hand. (See Sect. 7 for more details.) 

The present perspective, which is standard in condensed matter physics but has not yet been 
the subject of detailed mathematical analysis, can be summarized as follows: We will collect the 
important excitations about the various ground states into spin- wave modes. These modes form 
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the basis of an approximate low-temperature expansion which via the standard arguments yields 
the existence of several low-temperature states. 

Let us start by expressing all spins in terms of angular variables, i.e., Sr = {cos 9r, sin 6r). 
The (homogeneous) ground states are then just 6r = const = 6*. We define = Or — 6* so 
that, in the x direction, the interaction is given by ^(/3J)[cos(0* + ^r) — cos(0* + ^r+e^y? with 
analogous formulas in the y and z directions. Thus, to leading order in ^^'s, we have 



^4^Sf^ - ^Ssj' - ^^^\n{^r-irV.S (2.3) 



and, similarly. 



^ - S%f - ^ sin^(120° - r)(e. - irv..)\ 
^{SP - Sl%f « ^ sin2(120° + r)(e. - ^r+,f. 



(2.4) 



We will encode the 0*-dependence into effective coupling parameters: qi = sin^(^*), q2 = 
sm'^{9* — 120°) and = sin^(^* + 120°). Then the effective interaction for deviations about 
the ^*-state can be written as 

= ^ E E l-i^r - Cr+ej'. (2.5) 

a r 

Therefore, in some approximate sense, the partition function for deviations about the state where 
the spins are pointing in the direction 6* can be written as 

ZlA^*)^ |d^e-^^ef^(€) (2.6) 

where d^ denotes the product Lebesgue measure. 

As we will see, the integral is, as it stands, somewhat ill defined because the Hamiltonian pro- 
vides no decay for the zero Fourier mode of ^. However, it is recalled that for the above derivation 
to be meaningful, the S,/s had to be fairly small. So, one way out — which is what we will do in 
our proofs — is to restrict the integration measure in (2.6) only to (the Cartesian product of) small 
intervals centered at zero. Another way out, which leads to more transparent calculations, is to 
define the full object Z^A^*) the partition function constrained to configurations where, say, 
the average spin equals (cos 0*, sin 6*). (As we will see, inserting the appropriate (5-f unction on 
the right-hand side of (2.6) permits us to integrate the ^^'s over all real values.) In this language, 
the said constraint reads = 0, i.e., no "zero mode." For future reference, we denote the 

right-hand side of (2.6) with this constraint enforced by 



The reason for the prefactor will become clear momentarily. 



e-^^^^(^*). (2.7) 

V ^ J / 
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The translation-invariant structure of the effective Hamiltonian (and the constraint) prompts us 
to use the Fourier-transformed variables, 

r 

where A; takes values in = {27rL^-^(ni, n2, 713) : — L/2 < ni, 71,2,71-3 < L/2} which is known 
as the reciprocal volume (or the Brillouin zone). In terms of these variables, and the various other 
quantities defined, an appropriate spin-wave Hamiltonian can be constructed: 

*:eT*\{0} ^ a ^ 

where we have made explicit the absence of the contribution from the "zero mode." 
The calculation is now standard and we get 



n.h 



1/2 



(2.10) 



Thus, taking logs and letting L — ^ 00, we arrive at the limiting version of Fj^, 

This is the spin-wave free energy for fluctuations about the direction 0*. 

It is apparent that B* ^ F{6*) is invariant under the shift 9* ^ 9* + 60°. Far less obvious (but 
nevertheless true) is the fact that the absolute minima of F occur at 0* = 0°, 60°, 120°, . . . , 300°. 
Thus we must conclude that, when finite temperature effects are accounted for, six ground states 
are better off than any of the others. Sects. 4—5 will be devoted to a rigorous proof of this heuristic. 
A similar calculation allows us to estimate the spin-wave free energy for the inhomogeneous 
ground states and show that these are always less favorable than the homogeneous ones. 

Notwithstanding the appeal of the spin-wave heuristic, the above is just one step of the proof. 
In order to make use of spin- wave calculations, we resort to some (rather standard) contour esti- 
mates. Informally, we partition the "world" (by which we mean the torus) into blocks and mark 
those blocks where the spin configuration either features too much energy or has the character- 
istics of an environment without enough entropy. By adjusting the block scale we can make the 
penalty for marked blocks sufficient to carry out a Peierls argument. The principal tool for de- 
coupling the correlations between various boxes is provided by the chessboard estimates (which 
allow, via Cauchy-Schwarz type inequahties, to estimate the probabilities of various block events 
by their associated constrained partition functions). Explicit details are to be found in Sect. 6. 

Remark 3. It is noted that if the reader is willing to preaccept the forthcoming treatment as fact, 
an interesting feature concerning the surface tension is bound to arise. Indeed, let us imagine 
that the system is forced, e.g., via boundary conditions, to exhibit two favored states in the same 
vessel. The price for these circumstances will be the region — the interfacial region — where spins 
are bad. If /? » 1, the energetic form of "badness" can be ruled out a fortiori, but now we 
emphasize that the free energy difference between the most and least favored states is independent 
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of the temperature indicating that the cost of the interface will be temperature independent. Hence 
we anticipate that the stiffness (and also the correlation length) stays uniformly bounded away 
from zero and infinity as /3 — >^ oo. 

2.3 Discussion. 

The model under consideration exhibits infinitely many ground states, a problem which for math- 
ematical physics has surfaced but a few times in the past. When these situations arise, the finite- 
temperature fate of each ground state is typically decided by its capacity to harbor excitations. 
Here, the dominant excitations are exactly the spin waves from the last section — the spin-wave 
calculation shows that only a finite number from the initial continuum of ground states survive 
at positive temperatures. Unfortunately, an extra complication arises due to the inhomogeneous 
ground states discussed in Sect. 1.3. Here chessboard estimates allow us to bound the relevant 
spin-wave contribution by the spin-wave free energy against a periodic background. We remark 
that there are systems for which the spin-wave analysis featured herein may be performed with- 
out the complication of inhomogeneous ground states. One such example is the subject of the 
forthcoming paper [3]. 

As already noted, the "entropic-selection" mechanisms for long-range order are not new. In- 
deed, there have been some previous studies of the ANNNI models and other systems exhibiting 
infinite degeneracy of the ground state [7, 16,29]. However, the techniques involved in [7, 16,29] 
are based on the premise that there is a substantial gap in the energy spectrum which separates the 
excitations resolving the ground state degeneracy from the remaining ones. Due to the continu- 
ous nature of the spins, and the symmetry of the ground states, no such gap is of course present 
for the 120° -model. Instead, a decisive contribution to the entropic content comes from long 
wave-length excitations, i.e., the aforementioned spin waves. 

Another set of problems which are related to the present paper are the models with continuous 
spins studied in [17,49]. There the spins are a priori Gaussian random variables with covariance 
given by the inverse lattice Laplacian and with an on-site (anharmonic) potential. However, this 
potential is required to have only a finite number of nearly-quadratic minima (all of which have a 
uniformly positive curvature) which necessarily implies only a finite number of low temperature 
states. Notwithstanding, the work in [17,49] exemphfies situations where a ground state degener- 
acy is lifted by spin-wave-like excitations resulting in a reduced number of Gibbs states at positive 
temperatures. It is quite possible that the Pirogov-Sinai techniques used in [7, 16, 17,29,49] can 
after some work be adapted to our cases. However, at present the arguments via chessboard 
estimates seem considerably easier. 

As noted in Sect. 1 .2, the motivation to study these systems comes from the observed magnetic 
behavior of transition metal compounds. A complete understanding of these systems may there- 
fore require a full quantum-mechanical treatment. We expect a similar mechanism for ordering 
to be present also in the quantum-mechanical version of the 120° -model (as well as the orbital 
compass model). However, the only method of proof that seems promising in this context is the 
Pirogov-Sinai expansion of some sort. A general theory of these expansions for quantum systems 
exists, both for the situations with [15, 33] or without [6, 14] infinite degeneracy of the ground 
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state. But, as is the case for the classical systems, some fairly non-trivial generaUzations of the 
existing tools would probably be necessary. 

3. Proofs of main results 

In this section we will give the proof of our main theorem, subject to some technical results 
which will be proved later. In particular, in Sect. 3.1 we define the notion of a "bad" box and 
state without a proof the principal bound concerning the simultaneous occurrence of several bad 
boxes; see Theorem 3.1. This will be sufficient material for the proof of Theorem 2. 1 . The proof 
of Theorem 3.1 is the subject of Sects. 4-6; the actual proof comes in Sect. 6.3. 

3.1 Good and bad events. 

Here we will provide some mathematical foundations for our notions of the stable states and the 
contours that separate them. We will need three different scales — ^two of them spin-deviation 
scales and one a scale for the blocks which will be the setting of our various events. 

We will start with the fundamental spin-deviation scale which we denote by F. Here we are 
seeking a T (which is small) such that if all neighboring pairs of spins are within a distance F of 
each other, the harmonic approximation is "good" while if a neighboring pair violates this con- 
dition the energetic cost is drastic. On the basis of naive Taylor expansions — which is ultimately 
all we will do — it is clear that the latter is achieved if /3F^ » 1 and the former if /3F^ <^ 1. 
Thus, of course, we need P to be large and we can envision F to scale as any inverse power of (3 
between 1/3 and 1/2. 

The second deviation scale will be denoted by k and will serve to define sets of configurations 

which are effectively in one of the stable ground states. The third scale is the number B which 
will be used to define the spatial size of our block events. For fixed k, it appears that the only 
necessary requirement is that /3F^ » log B which will always hold eventually. Unfortunately, 
there is some spurious interplay between the parameters B, k and F which could, in principle, be 
removed in a more refined analysis. But, for this work, we will keep the "smallness" of k in the 
realm of the existential and require B to get large, but only very slowly, as (3 goes to infinity. 

In order to make our main technique, the chessboard estimates, available we have to con- 
fine ourselves to systems with periodic boundary conditions. Let thus T^^ denote the three- 
dimensional torus of scale L, i.e., = 1? jiUL)^ . In general, we will be dealing with certain 
events taking place in blocks of a specific scale B and we will be using the chessboard estimates 
to bound probabilities of these events. These blocks will be translates of the block As C 
which we define as the cube of {B -|- 1)^ sites with the "lowest left-most" site at the origin. It will 
be convenient (although presumably not strictly necessary) to assume that the Unear scale of our 
finite-volume system, L, when divided by B results in a power of two. 

Now we are ready to state the definition of a "good" block: 

Definition 1 Let B denote a positive integer and let /t > and F > be sufficiently small. We 
will say that the spin configuration in the block A b (or the block itself) is good if the following 
two conditions are met: 
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(a) 



For each a G {1, 2, 3} and any neighboring pair r and r + e^ in A^, 




(3.1) 



(b) All spins in A^ make an angle which is less than 2k from one of the preferred six directions. 
Exphcitly, if Sr = (cos 0^, sin 9^) then, for some r = 1, . . . , 6, we have \dr — < 2k for 
all r e A^. Here, of course. Or is only determined modulo 2it. 

We denote by ^ = Gb,k,t the event that the block A^ is good. The complementary event, 
marking the situation when the block is bad, will be denoted by inevitable B. Our goal will be 
to bound various probabilities involving bad events. The main tool for these bounds will be the 
chessboard estimates whose basic setup and principal result we will now describe. 

Let Fl,i3 and (— denote the (Gibbs) probability measure, respectively, the corresponding 
expectation according to the Hamiltonian (1.1) at inverse temperature /3 on T^. Let t denote a 
vector with integer coefficients identified modulo L / B — in formal notation, t G T^,/ ^ — and let B 
be an event discussed above. Then we let 'dt{B) denote the event B translated by the vector Bt. 
(For general events A defined on the configurations in A^ we will need an enhanced definition 
of ??t(.4); cf the definition prior to Theorem 6.2.) Note that if '&t{B) and {B) are "neighboring" 
translates of B, then these two events both depend on the spin configuration on the shared face of 
the corresponding translates of A^. 

The principal result of this section, which is the starting point for all subsequent results of this 
work, is the following theorem: 

Theorem 3.1 Consider the 120° -model as defined by (1.10). For each sufficiently small k > 
and each rj £ (0, 1) there exist Lq G (0, oo) and (3o G (0, oo) and, for any any (3 > f3o, there 
exist numbers T G (0, 1) and B G (0, oo) such that the following holds: If B is the event — 
defined using k, T and B — that the configuration in Kb is bad and ti, . . . ,tm are distinct vectors 
from T^/s, then for any L > Lq, 



This result provides a way to estimate the probability of simultaneous occurrence of several 
bad events. The non-trivial part of the proof of Theorem 3.1 boils down to the spin-wave calcu- 
lations outlined in Sect. 2.2. The rigorous version of these calculations requires some substantive 
estimations and the actual proof is therefore deferred to Sect. 6.3. 

3.2 Proof of long-range order. 

Now we are ready to prove our main theorems. We note that there are six disjoint ways to exhibit 
a good block for the 120° -model, each corresponding to one of the vectors w,-. We will denote 
the corresponding events by Qr, with r = 1, 2, . . . , 6. Exphcitly, 



PL,/3K(i3)n---n^t^(i3))<r?- 



(3.2) 



gr = gn{S: Sr- 



Wr > cos(2k), r G A^}. 



(3.3) 



The core of the proof is the following (almost direct) consequence of Theorem 3.1: 
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Lemma 3.2 Consider the 120° -model onT^ and suppose that k <^ 1. There exists a function 
h : [0, 1) — > [0, oo) satisfying h{rj) as rj I 0, such that for each sufficiently small rj > and 
each j3, T and B as allowed by Theorem 3. 1 the following is true: For any t\,t2 G '^l/b cind any 
type of goodness t, we have 

nA^tAGr) n ^t,{Sr)) < h{v), (3.4) 

provided L > Lq, where Lq = Lq{k, rj) is as in Theorem 3.1. 

Proof. Noting that Q^- nQ^ = 0, let us assume that ti 7^ t2- Now, for the intersection {Gr) H 
't?i2 (St) to occur, either the block at Bt2 is bad, which has probability at most t], or it is good but 
not of the type r. We claim that, in the latter case, there must be a "surface" consisting of bad 
blocks which separates the block at Bti from that in Bt2. Indeed, let S G -dti (Gt) and consider 
the connected component, ^, of good blocks in Tx,/^ containing the block at Bti. We claim that 
the type of goodness is constant throughout i.e., it is of type r. To see this, suppose that a 
block in ^ has the type of goodness which is distinct from r. By the fact that is connected, 
there must exist a pair of neighboring blocks with distinct types of goodness. But neighboring 
blocks share the sites on their separating face and (since k -C 1) the spins on this face cannot 
simultaneously be in the 2K-neighborhood of two Wt-'s — that is, not without the spins busting 
apart. Hence, on 1?^^ (Gt) n ■dt2 (Qr)' ^^e block at Bt2 is not part of ^ and we have it separated 
from Bti by a (*-connected) "surface" of bad blocks. 

To estimate the probability of such a "surface" we will use Theorem 3.1: The probability that 
a "surface" involving altogether m given bad blocks occurs is bounded by rj"^. The rest of the 
proof parallels the standard Peierls argument which hinges upon the fact that the number A^^, of *- 
connected "surfaces" comprising m blocks and containing a given block grows only exponentially 
with m, i.e., < for some c G (1, 00). To count the number of ways how to choose the 
particular block in the "surface," we have to be a bit cautious about the toroidal geometry: If 
m < L/B, then the "surface" encloses either the block at Bti or that at Bt2 on all sides and 
there are at most 2m ways to choose one particular block. On the other hand, if m > L/B, 
then the surface can be topologically non-trivial but, since Ti/j^ is a finite graph, the number of 
choices of one particular block is at most (L/B)^ < m^. This shows that (3.4) holds with 

h{r]) = + 2 ^ m^{cr])'^, (3.5) 

m>6 

uniformly in L > Lq. Clearly, h{ri) ^ as 77 | 0. □ 

Now we are ready to prove the existence of long-range order in 120° -model, subject to the 
validity of Theorem 3.1: 

Proof of Theorem 2.1. Let 77 > and let /3o and Lq be as in Theorem 3.1. Fix a P > Po 
and choose B and T accordingly. For finite L > Lq, it follows by (3.2) that, with probability 
exceeding I — r], any given block is in a good state. Since the distinct types of goodness are 
disjoint and related by symmetry, we have 

^LAMOr))>kl-v) (3.6) 
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for any t G T^/^ and any r = 1, 2, . . . , 6. Next, we may condition the block farthest from the 
origin (i.e., the one at the "back" of the torus) to be of a particular type of goodness, say Qr- The 
resulting measure still satisfies the DLR-condition in any subset of the torus not intersecting the 
far-away block. Passing to the thermodynamic Umit along some sequence of L's, we arrive at an 
infinite- volume Gibbs state for the interaction (1.1). 

Clearly, by (3.6) and Lemma 3.2, we have the uniform bound 

rLA^tAG'r)K{Gr))<Q^. (3.7) 

Hence, if <C 1, we have constructed six infinite-volume Gibbs states in the 120°-model which 
are distinguished by the statistical properties of any individual spin. In particular, the bound (2.2) 
holds with e{P) directly related to rj, h{ri) and n. Of course, it is not automatically the case 
that the resulting states are translation-invariant; however, this is easily handled by considering a 
translation average of the abovementioned and noting that the "distinctness" of the states via the 
single spin observables is preserved by this averaging. □ 



4. Spin- WAVE analysis 

This section provides rigorous justification for the heuristic spin-wave calculations from Sect. 2.2. 
Beyond the fact that these calculations settle the pertinent questions concerning long-range order 
at the non-rigorous level, such results, as refined here, serve as the cornerstone for the proof of 
Theorem 3.1. The principal results of this section are Theorems 4.1 and 4.5. 

4.1 Homogeneous ground states. 

Our goal is to evaluate the free energy of the spin configurations where all spins are more or less 
aligned with a given vector on the unit circle. Let us represent all of the spins S = (Sr) by their 
corresponding angle variables 6 = {Or) — vis-a-vis the usual Sr = {cos 6r, sin 6r) — and let 6* 
denote the particular direction towards which we wish the spins to align. Let xa,l(^) be the 
indicator of the event that j^^ — < A, with the difference Or — 0* interpreted modulo 27r, holds 
for all r G T^. Here A is closely related to the quantity T from Sect. 3.1. 
In this representation we define the constrained free energy by the formula 

where J^l{0) denotes the torus Hamiltonian expressed in terms of the angle variables and 
where the first term on the right-hand side has been added for later convenience. Our goal is 
to show that, under specific conditions, {0*) can be well approximated by the function F 
defined in (2.1 1). (As is easy to check, the integral in (2. 11) converges for all 0*.) 

Recall the abbreviations gi = sin^(0*), q2 = sin^(0* - ^) and = sin^(6'* + ^) from 
Sect. 2.2. The precise statement concerning the above approximation is as follows: 
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Theorem 4.1 For each e > there exists a number 5 = 5{e) > such that ifPJ and A obey 

{PJ)A^ > 1/S and {pJ)A^ < 5, (4.2) 

then 

lim sup I F^^^ {e*) - F{e*) I < e (4.3) 

for all e* e [0,27r). 

As the first step of the proof, we will pass to the harmonic approximation of the Hamiltonian, 
which is given by 

MO) = ^ E E ^«(^' - ^^+6.)'- (4.4) 
The next lemma provides an estimate of the error in this approximation. 

Lemma 4.2 There exists a constant c\ G (0, oo) such that for all /3 G (0, oo), all A G (0, 1), 
all L > 1 and all 0* € [0, 27r) the following holds: IfxA,L{0) = 1, then 

\pj^L{e) - J^L{e)\ < ci{pj)A^L\ (4.5) 

Proof. Let us first consider the nearest-neighbor bond (r,r + ei) and note that Sr^^ = cos 9r. 
Since \6r — 0*\ < A, Taylor's Theorem gives us the bound 

\S^'^ - S^%^ + sin(r)(e, - ^,+ei)| < (4.6) 

But {Or - 0;.+eJ < 2A and thus (5^^^ - 5"^+^^^ and qiiOr - 9r+g^)^ differ by less than a 
numerical constant times A'^. The situation in the other directions is similar, one just has to 
note that sP = cos{e - ^) and Sr^^ = cos{9 + Adding up the contribution of all three 
components, multiplying by (3 J and summing over r G T^, the result directly follows. □ 

Having converted the Boltzmann weight q-P-^l^S) ^^i^^ Gaussian weight e~'^^^^^ in (4.1), 
our next task is to estimate the effect of the indicator xa,l- Let 

Then we have: 

Lemma 4.3 For all P G (0, oo), all A G (0, 1) and all 9* G [0, 27r), 

lim sup g'^ < -F(r). (4.8) 

Proof. We will use the exponential Chebyshev inequality. Let A > 0. Then the indicator xa,l is 
bounded via 

Xa,l(0) < e^^(/^^)^^^^ exp{-iA(/3J) ^ (^^ " n'}- (4-9) 
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Plugging the right-hand side into (4.7) instead of xa,l. we get a Gaussian integral with L^- 
dimensional covariance matrix C = {|3J)~^{X^ + D)~^, where 1 is the unit matrix and D is a 
generalized Laplacian implicitly defined by (4.4). Integrating out the variables 6 and invoking 
Fourier transform to diagonahze C, we get 

< IxmA^ - li, ^ log{A + (4.10) 

JteT* a 

where denotes the reciprocal volume (or the Brillouin zone). Passing to the limit L — > oo, we 
find out that the left-hand side of (4.8) is bounded by ^X{/3J)A^ - F{e*, A), where 



F{9 



But the integrand is a monotone function of A, and so the Monotone Convergence Theorem 
guarantees that F{9*, A) J, F{9*) as A J, 0. Thence the result follows by taking A to zero. □ 

Let F{9* , A) be the quantity defined in (4.11). The lower bound is then as follows: 

Lemma 4.4 For all j3 G (0, oo), all A G (0, 1), all 9* G [0, 2n), and all X > satisfying 
{(3J)A'^X > 1, 



iim ini 



........ ^3 >-f(«.,A)+log(l-^i). (4.12) 

Proof. Let A > and consider the Gaussian measure Pa given by 

where Q L,e* is the corresponding normalization factor, which modulo the "log" and the factor 
equals the final term on the right-hand side of (4.10). Let us use Ex to denote the corresponding 
expectation. Then 

QfT = QLfi*mXA,L)- (4.14) 

Now XA,L is simply the product of indicators of the type 1 {|er-6i*|<A}- We claim that 

mXA,L)=Ex[ll ^{ler-e*\<A}) > n M\^r-9*\<A). (4.15) 

This follows from the fact that the moduli of these sorts of Gaussian fields are FKG-positively 
correlated, see e.g. [5]. It can also be established on the basis of the "esoteric" version of reflection 
positivity (using reflections between sites), which is described at the beginning of Sect. 6.1. The 
estimate thus boils down to a lower bound on the probability of 1^?^ — < A. 

Now, let us note that the Fourier components 9k of the fields 9r — 9* have Ex{9k) = and, 
for k' 7^ ±k, the random variables 9k and 9k' are independent with 
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Thus, invoking the Chebyshev inequality the complementary probability is bounded by 

F,{\er -e\>A)< ^ 2^ < (4.17) 

Combining (4.14-4.15) with (4.17), invoking the explicit expression for (5l,6»* and passing to the 
Umit L — >^ oo, the desired bound is proved. □ 

Now we are ready to prove the error bound in (4.3): 

Proof of Theorem 4.1. By the Monotone Convergence Theorem we have that F{9*, A) J. F{9*) 
as A 1 0. Moreover, the continuity of 0* i-^ F{0*) and the fact that the unit circle in R2 is 
compact imply that this convergence is actually uniform in 9*. Hence, for each e > 0, there exists 
a number A > such that 

\F{9\\)-F{e*)\<^- (4.18) 

for 9* G [0, 27r). Let ci be the constant from Lemma 4.2 and choose 5 such that ci5 < e/3. 
Suppose also that S < X and 

log(l4)>-l (4.19) 
Fix an angle 9* G [0, 27r). Lemma 4.2 along with our choice of S imply 



lim sup 



< |- (4.20) 



L3 ^ ^ L,f3 

On the other hand. Lemmas 4.3^.4, the choice of A in (4.18) and our choice of S ensure that 



lim sup 

L—too 



< (4.21) 



L3 

Combining these two estimates, the bound (4.3) is proved. □ 
4.2 Stratified ground states. 

As mentioned previously, constant configurations are only the overture for the set of all possible 
ground states. As a consequence, the knowledge of the spin- wave free energy about homogeneous 
background configurations is not sufficient for the proofs of our main results. Fortunately, as we 
shall see in Sect. 6, the chessboard estimates allow us to reduce the (potentially quite large) 
number of remaining cases to configurations which are translation-invariant in two directions and 
alternating in the third direction. 

To avoid parity problems, throughout this section we will assume that L is an even integer. Fix 
an index a € {1, 2, 3}, pick a direction 9* G [0, 27r) and let 9* denote the reflection of 9* through 
the a-th of the vector a, b or c. Consider again the angle variables 9r and let xa,l be the indicator 
that \9r — 9*\ < A f or r G Tl with an even a-th component while for r with an odd a-component 
we require that \9r — 9*\ < A. Let 

Fi^'-'in = log - ^ log / e-'^'^^'^XAAO) n d^^- (4.22) 



Foc{0*) = -,1 7:^logdetn„(A:). (4.25) 
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The quantity F^p"'\6*) represents the spin-wave free energy for (period-two) stratified states 
perpendicular to direction a and spins alternating between directions 0* and 9*. 

As before, our goal is to approximate {0*) by an appropriate momentum-space integral. 

For a G {1, 2, 3}, let us abbreviate 

Ea = Ea{k) = \l - e'''"\^ and E% = E'^{k) = \l + e'^"\^ , (4.23) 

where is the a-th component of the vector k G [— tt, tt]^, and recall, once again, the meaning 
of the quantities (cf Sect. 2.2). We will define three 2 x 2-matrices Yla{k), a = 1,2, 3. First 
let a = 1 and abbreviate g+ = ^{q2 + Qs) and g_ = ^{Q2 — Qs)- Then 

V Q-{E2-Es) q,E*,+q+{E2+Es)J- ^^"^^^ 

The quantities 112 and 11,3 are defined by cycHcally permuting the roles of Ei , E2 and E^ and 
similarly for the g^'s. (In the physically relevant quantities, g_ will appear only in terms of its 
square, so the order used for the definition of this quantity is for all intents and purposes arbitrary.) 
Then we define a function F^ assigning to each 6* G [0, 2tt) and each a G {1, 2, 3} the value 

1 r dk 

The fact that now we have a quarter in front of the integral comes from the fact that the determi- 
nant actually represents the combined contribution of two A: -modes. 

The main result of this section concerning the 120°-model is now as follows: 
Theorem 4.5 For each e > there exists a number 5 = 5{e) > such that if PJ and A obey 

{PJ)A^ < 6, (4.26) 

then 

lim inf F^J);"^ {9*) >F^{e*)-e (4.27) 
for all 9* G [0, 27r) and all a e {1,2, 3}. 

As we have seen in the previous sections, the first step is to pass to the quadratic approximation 
of the torus Hamiltonian. Let = {9r)'bea configuration of angle variables. Then we define 

^lAO) = ^ E E '^feV^ - ^^+e,)^. (4.28) 

reTi 7=1,2,3 

(a) 

Here q^J is the usual q^ if the cc-th component of r is even while for the complementary r we 
have to interchange the roles of the two qy with 7' 7^ a. (In particular, q^} = q^ for all r.) 
Our first item of concern is the error caused by this approximation: 

Lemma 4.6 There exists a constant C2 G (0, 00) such that for all (3 G (0, 00), all A G (0, 1), 
all L > 1 and all 9* G [0, 27r) the following holds: If xa,l{0) = 1 and if = {9r) is the 
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configuration obtained by reflecting 9r through the a-th of the vectors a, b or c for r with an 
odd a-component, then 

\pJfL{e) - ^lAH < C2{f3J)A^L\ (4.29) 

Proof. Once we have accounted for the inhomogeneity of the setup, the proof is essentially 
identical to that of Lemma 4.2. Without much loss of generality, let us focus on the case a = 1. 
Let be such that xa,l(^) = 1 and let 6 be as described. 

We will concentrate on the interaction of two spins in the y-coordinate direction. If the a- 
component of r is even, then the expansion around 9r = Or ^ 0* gives that cos{6r — ^) is well 
approximated by cos{9* — ^) — sin(0* — ^){9r — 9*). Accounting better for the errors we 
thus get 

1-5?^ - <^e. + - f )(^. - er+i,)\ < (4.30) 

On the other hand, for r with an odd a-th component we have —9r = 9r^9* which means that 
SP = cos{-9r - ^) = cos{9r + ^) and thus 

1^/'^ - Sl%^ - sin(r + ^){9r - < A^. (4.31) 

After plugging into J^l (6), the 9r in the even r planes are coupled by q2 while in the odd r planes 
they are coupled by qs, in accord with (4.28). 

A completely analogous argument handles the case of two sites in the z-coordinate direction. 
In the x-direction the reflection has no effect because the minus sign from sin disappears 
after we take the square. The ensuing errors are estimated exactly as in the proof of Lemma 4.2. 

□ 

Since xa,l{0) = XA,L{d), where 6 and 6 are related as in Lemma 4.6, a simple change of 
variables shows that the proper analogue of the quantity from (4.7) for the present setup is 

= (1^)'^'' / --^~^-'''xaAO) n ^^r. (4.32) 

Note that here the inhomogeneity of the domain of integration in (4.22) has now been moved to 
the Gaussian weight. Next we apply: 

Lemma 4.7 For all (5 G (0, oo), all A G (0, 1), all 9* G [0, 27r) and all a G {1, 2, 3}, 

limsup ^3'^'^ < -Fai9*). (4.33) 

Proof. Fix a and let A > 0. The proof again commences by invoking the exponential Chebyshev 
inequality in the form of (4.9). The resulting -dimensional Gaussian integral has covariance 
matrix Ca = (/3J)(A1 + D^"'^), where ^f3JD^"'> is the matrix corresponding to the quadratic 
form (4.28) in the variables 9r. The difference compared to Lemma 4.3 is that now is no 
longer translation-invariant in the a-th direction, but only periodic with period two. As a result. 
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the k and k + TrCo, modes will mix together and the Fourier transform of will result in 2 x 2- 
block-diagonal matrix. The blocks are parametrized by pairs of momenta {k,k + ttSq). (Note 
that, since L is even, k + ttSq G whenever k G 

A calculation — which is best performed by taking the Fourier transform of ,-^L,a — reveals that 
the block corresponding to the pair {k,k + ire-a) is exactly Ila{k). Hence we get 

°^ J'^'^ < ^A(/?J)A2 logdet(A1 + U^{k)), (4.34) 

where 1 is the 2 x 2-unit matrix and where the usual factor 1 /2 in front of the sum is replaced 
by a 1 /4 to account for the fact that k and k + ttCq are treated as independent entities in the sum. 
(We are using the A: A: + ttCq symmetry of the determinant.) Passing to limits L ^ oo and 
A i 0, the bound (4.33) is proved. □ 

Proof of Theorem 4.5. By Lemma 4.6 and the definition of F^^f'' (9* ) we know that 

Fl^/)(r) > -i^^l^M _ c,{(3J)A\ (4.35) 
Hence, if S is such that C2S < e, (4.27) follows by taking L ^ 00 and invoking Lemma 4.7. □ 



5. Spin- WAVE free energy minima 

The purpose of this section is to show that the spin-wave free energy F{9*), which emerges from 
the analysis in Sect. 4.1, is minimized in the "directions" as stated in Theorems 2.1. Similarly, 
we will also show that the free energy Fa {9*) corresponding to the inhomogeneous ground states 
is always strictly larger than its homogeneous counterpart F{9*), unless 9* is "aligned" (or "an- 
tialigned") with the a-th of the vectors a, b or c. These findings constitute the essential ingredients 
for the proof of Theorem 3.1 in Sect. 6.3. The principal estimates are based on Jensen's inequality 
combined with a non-trivial bit of "function analysis." 

5.1 Homogeneous ground states. 

Our task is to identify the minima of the function 1-^ F{9*) defined in (2.11). However, 
noting that the product structure of the measure dk / (27r)^ makes the random variables 1 1 — e''^" | ^ 
independent, we might as well analyze an entire class of functions of this type. 

Let X be a random variable taking values in [—1, 1] and, for any triple of numbers (a, b, c), 
define the function 

F{a, b, c) = E(log(a2(l - X^) + b''{l - X2) + c\l - X3))) , (5.1) 

where Xi,X2 and X3 are independent copies of X. Suppose in addition that the distribution n 
of X has the following properties 

(1) n has a density with respect to the Lebesgue measure, /x(dx) = f{x)dx. 

(2) /(x) = /(-x)forallxe [-1,1]. 
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(3) f{x) is Strictly increasing on [0, 1]. 
Then we have the following general result: 

Theorem 5.1 Let a,b,c i-^ F(a, b, c) be as above with X satisfying the properties (1-3). Then 
for any nonzero x G M and any o, 6, c satisfying 

a + 6 + c = and + 6^ + = 2x^, (5.2) 

we have 

F{a, b, c) > F(0, X, -x). (5.3) 
Moreover, the inequality is strict whenever a,b,c^ and F{0, x, — x) > — oo. 

The particular case of the 120° -model can now be easily extracted: 
Corollary 5.2 The function 9* i-^ F{9*) achieves its minimum only at the points 

^* = |r, T = l,2,...,6. (5.4) 

Proof. We just have to identify the quantities a, b, c and the random variable X in the case of 
the 120° -model. First, since |1 — e'^"p = 2(1 — cos we let X be the random variable 
distributed as cos/c in measure dfc/(27r) on [— 7r,7r]. A simple calculation shows that X has a 
density f{x) = (1 — x^)~^/^/(27r) with respect to the Lebesgue measure on [—1,1], which 
manifestly satisfies the requirements (2-3) above. 

Now, setting a = V2sm{e*), b = v^sin(6'* - ^) and c = V2sm{e* + ^), we have 
that F{a, b, c) = 2F{9*). Moreover, a trivial calculation shows that a + b + c = 0, while 

+ b'^ + c'^ = S and (5.2) thus holds with x^ = 3/2. As a consequence, 9* hh- F{9*) is 
minimized only by 9* such that one of the numbers a, b, c vanishes. This is easily checked to give 
just the values in (5.4) □ 

The rest of this section is devoted to the proof of Theorem 5.1. The proof is based on two 
observations: First, a lemma due to Pearce [40]: 

Lemma 5.3 Let X be a random variable on [—1,1] satisfying properties (1-3) above. For 
each A G M, let {—)\ denote the expectation with respect to the probability measure u\{dx) = 
Nxf{x)e^^dx, where f is the probability density of X and Nx is an appropriate normalization 
constant. Then the function A ^ {X)\ is strictly concave on [0, oo). 

Proof. See [40]. The conditions (1-3) represent one of the general situations in which one can 
prove the GHS inequality in lattice spin systems, see [44, Theorem n.l3.5(ui)]. □ 

The second observation is that the constraints (5.2) reproduce themselves, rather unexpectedly, 
at the level of quartic polynomials in a, b and c: 

Lemma 5.4 Let a, b, c be numbers satisfying (5.2). Then 



a^ + h^ + c'^ = 2x\ 



(5.5) 
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Proof. Since o = — (6 + c), eliminating a from the second constraint in (5.2) results in 

+ + bc = x^. (5.6) 

Squaring we get 

b^ + c^ + b'^(? + 2b'^c^ + 26c(62 + (?) = x^ (5.7) 
which can be recast into the form 

26^ + 2c* + eft^c^ + 463c + 4hc^ = 2x^. (5.8) 
Splitting off the term b* + c^, the rest of the left hand side is clearly {b + c)^ = a*^. □ 
With these lemmas in the hand, the proof of Theorem 5.1 is relatively straightforward: 

Proof of Theorem 5.1. Since we can scale a, b and c by any constant at the cost of changing 
F(a, 6, c) only by an additive x-dependent factor, let us suppose without loss of generality that 
X = 1/ \/2. Moreover, if one of a, 6, c is zero, say a = 0, then h = — c = ibx and (5.3) is trivial. 
Hence, we only need to focus on the situations when a, 6, c / 0. 

The first step of the proof is to convert the expectation of the logarithm into the expectation of 
an exponential. This can be done for instance by invoking the identity 

r-l r-oo -s 

-log(l-x)=/ di / ds— (e"*^-l), x < 1, (5.9) 
Jo Jo t 

where the double integral on the right-hand side is well defined because everything is positive. 
Let us now plug in a?Xi + 6^X2 + c^X^ for x and take expectation with respect to Xi, X2 
and X3. Applying Fubini's theorem (and the fact that, almost surely, a^Xi + 6^X2 + c^X^ < 1), 
the result is 

f-l f-OO -s 

F{a,b,c)= dt ds {l - G{st;a,b,c)) , (5.10) 

Jo Jo t 

where 

G(A; a, b, c) = E{e^^''^^^+^^^^+^^^^^). (5.11) 

We will show that, whenever a,b,c ^ 0, we have G{\; a, b, c) < G(A; 0, l/\/2, — 1/\/2) for all 
A > 0, from which (5.3) and the ensuing conclusion directly follow. 
Consider the function A 1— > R{X) defined by 

i?(A) = log GiX;ab,c) 

' ^G(A; 0,1/^2,-1/^/2) 

Our goal is to prove that R{X) < whenever A > 0. First we note that R{0) = so it suffices to 
show that R'{\) < for all A > 0. Invoking the independence of Xi, X2 and X^, we have 

R'{X) = a^{X)^,2 + b^{X)^,2 + ?{X)^,2 - {X)^/2, (5.13) 

where we adhere to the notation from Lemma 5.3. Now, by Lemma 5.3 and our assumptions on 
the random variable X, the function A <—>■ {X)\ is strictly concave. Since + b^ + (? = 1 and, 
as guaranteed by Lemma 5.4, also a* + b* + = 1/2, the bound R'{X) < for A > is a direct 
consequence of Jensen's inequality. □ 
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Remark 4. The previous proof has one (arguably) unnatural feature; namely, the conversion 
"from logs to exponentials" via the identity (5.9). It would of some interest (at least for the 
authors) to see whether a more direct argument can be constructed. 

5.2 Stratified ground states. 

Having identified the absolute minima of the spin-wave free energies for homogeneous back- 
ground configurations we turn our attention to the free energies corresponding to inhomogeneous 
ground states. Specifically, we will show that (truly) stratified states have always worse free 
energy than the corresponding homogeneous ones. 

Let F{9*) denote the spin-wave free-energy from Sect. 4.1 and let Fa{9*) be the corresponding 
quantity for the stratified states as defined in Sect. 4.2. Then we have: 

Theorem 5.5 For each k > there exists a constant C3 = c^{k) > such that if a G {1, 2, 3} 
and if the angle between 6* G [0, 2tt) and the a-th of the vectors a, b and c is in {k, tt — k), then 

Fa{0*)>F{e*) + C3. (5.14) 

Proof Recall the notations Fa = |1 — e'*^"|^, F* = |1 + e'*^"p and qi = sin^(6'*), q2 = 
sin'^ie - ^) and 53 = sin2(6' -h ^) and the definition of Ua{k) in (4.25). We will write 
det Haik) as a convex combination of two terms each of which produces the same free energy. 
Without loss of generahty, let us assume that a = 1. We claim that for all k G [— 7r,7r]^, the 
quantity Hi (A:) admits the decomposition 

det Ui{k) = ^ {qiFi + 52^2 + 93^3) {qiEt + 93^2 + q2E3) 

+ ^ {qiEi + q3E2 + g2£^3) (gi^? + q2E2 + qsEs) ■ (5.15) 

To prove this let us abbreviate q± = 5(^2 i §3) and E± = £^2 =b £'3. Focusing on the first term 
on the right-hand side, we write 

qiEi + 92-^2 + 93-^3 = qiFi + q+F+ + q-E_, 

qiE* + 93^2 + 92^3 = qiEt + q+E+ - q^E^, ^^'^^^ 

and similarly for the other two terms. Multiplying these two lines tells us that the first term on 
the right-hand side of (5.15) equals a half of 

(gi^l + q+F+){q^F* + q+F+) - q^Fl + qiq-{F\ - F^)F_. (5.17) 

The sole effect of the second product on the right-hand side of (5.15) is to cancel the very last 
term of (5.17) — ^note that the sign of g_ changes when q2 and ^3 are interchanged. Now the first 
two terms in (5.17) is exactly the determinant of Hi (A:). Hence (5.15) follows. 

If we plug in any of the four linear factors in Fi, F2.. E3 on the right-hand side of (5.15) into 
the logarithm in (4.25), integrate and apply the symmetries of the measure dk, the result will be 
lF{e*). Suppose now that 9* / 0°, 180° and note that this implies that q2 ^ qs- Then 92 £^2 + 
qsEs ^ qsE2 + 92 ^^3 on a set of positive measure dk. Hence, the two terms on the right-hand 
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side of (5.15) are not equal almost surely which by the strict concavity of the logarithm and 
Jensen's inequality implies that Fi{6*) > F{6*). But both functions are continuous in 9*, and so 
Fi{9*) — F{9*) is uniformly positive on any compact subset of the unit circle not containing 0° 
and 180°. The existence of a desired C3 follows. □ 



6. Probabilities of bad events 

Our goal here is to prove the estimate in Theorem 3.1 concerning the probability of a simultaneous 
occurrence of several bad events. While some of the details may still appear to be rather intricate, 
the principal input into the forthcoming argument has already been established in Sects. 4-5. 

6.1 Reflection positivity and chessboard estimates. 

In this section we will glean from the classic theory of reflection positivity those items that are 
needed at hand. Recall our notation g for the Gibbs probability measure on defined by the 
Hamiltonian (1.1) at inverse temperature 13. Reflection positivity is a property of the measure ^l,i3 
with respect to reflections of the torus which are defined as follows: Suppose that L is even and 
let us split the torus symmetrically into the "left" and "right" parts, and TJ, such that the 
two reflection-symmetric halves either share two planes of sites (reflections "through sites") or 
not (reflections "through bonds"). Let P be the formal notation for the "plane of reflection" 
and let ^p, resp., denote the cr-algebra of events that depend only on the portion of the 
configuration in T^, resp., T£. Introduce the reflection operator ■i^p on configurations in T^, 
which induces a corresponding map d p : ^p — > . Then we have: 

Lemma 6.1 (Reflection positivity) Consider the plane P, the cr-algebra and the mea- 
sure Pi ,3 as specified above. Let 1El,/3 denote the expectation with respect to Pl,/3- Then the 
following holds for all bounded ,^p -measurable random variables X and Y: 

^L,p{X^p{Y)) =^L,p{Y^p{X)) (6.1) 

and 

El,/3(^^?p(^)) > 0. (6.2) 
Here tI}p{Y) denotes the random variable Y o i^p, and similarly for '&p{X). 

Proof. This is the standard reflection positivity proved in [22-24], which for reflections "through 
sites" follows simply by the fact that the interaction is exclusively via nearest neighbors, while 
for reflections "through bonds" it follows from this and the fact that the coupUng is both quadratic 
and attractive. □ 

Remark 5. We remark that in the present work we use only the more robust version of reflection 
positivity — ^poor man's RP — which only requires nearest-neighbor interactions. (An exception 
to this "rule" is perhaps the argument leading to (4.15); but there we also offer an alternative 
approach via [5].) 
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Our use of reflection positivity will come through the so called chessboard estimates. To 
motivate the forthcoming definitions, let us briefly recall the principal idea. Using the expression 
on the left-hand side of (6.2), one can define an inner product on the -measurable functions, 
which then satisfies the Cauchy-Schwarz inequahty 

PL,^(^n^p(^'))' < PL,/3(^n^?p(^))P^,^(^'n7?p(^')), (6-3) 

for any A,A'^ ^p. The interpretation of this inequality is that two given events, one on the 
"left" and the other on the "right" of T^,, can be separated within the expectation at the cost of 
reflecting both of them through P. Iterating this bound further one can eventually disseminate 
each event all over the torus. The resulting quantity is often amenable to further analysis. 

To state the chessboard estimates formally, let us consider a rectangular box F c of 
dimensions ai x 02 x 03, where the a^'s are positive integers. For simpUcity, here and throughout 
this work, we assume that all of the ai are related to L by powers of two, i.e., ai = 2~™'L for 
some integers m^. Consider the tihng of the (continuous) torus with dimensions L x L x L by 
translates of V. We will parametrize these translates by vectors teT where T is the (discrete) 
torus with dimensions L/ai x L/a2 x L/a^. 

Let A be an event which depends only on configurations in F n T^. First we note that the 
event A can be reflected (multiply) through the various midplanes of V , leading to seven new 
ostensibly different versions of the event A. [Labeling the resulting events by cr = (cji, 02, 0-3) G 
{0, 1}'^, where 0-^ = 1 denotes whether the reflection in the a-th direction is implemented, we 
thus have altogether eight events: one ^000 = -4, three order- 1 reflections Aiqq, Aqiq and ^001 
through the midplanes of V orthogonal to x, y and z lattice directions, respectively, three order-2 
reflections .Aiio,^ioi,-4oii and one order-3 reflection ^m.] Nowift G T, let us define 
the appropriate notion of "translation by t," as follows: For t's with all even coordinates, this is 
simply the usual translation by t. For t's with some odd coordinates, we select from the other 
seven versions of A the one with reflections corresponding to all the odd coordinates of t; the 
event 'dt{A) is then the translation by t of that version of A. 

Let Zl^p{A) denote the constrained partition function defined by 

ZI^,p{A) = ZL,p{X{^,,^A))^^^, (6.4) 

teT 

where Zl^i^ is the usual partition function on Tx, and denotes the indicator function of 

event "dtiA). We are now ready for: 

Theorem 6.2 (Chessboard estimate) Let the events Ai,..., Am and the partition functions Zj^^jj 
and ZL^jj(^Ak) be as above and let ti, . . . ,tm be distinct vectors of the type described. Then 

PM(n*..(A))<n(^f^)"*. (6.5) 

k=l k=l ^'^ 

where \T\ is the volume of the factor torus T. 
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Proof. This is the standard chessboard estimate proved in [22-24], see also [43]. These estimates 
follow, in general, whenever the interaction is reflection positive — here using both reflections 
through sites and bonds depending on whether the corresponding coordinate of the side of V is 
integer of half-integer. □ 

Remark 6. We observe that due to the self-imposed evenness constraint on the dimensions of the 
tori, the objects ZL^f){A) and ZL^p{'dt{A)) are identical for all t. This will reduce the need for 
various provisos in the future derivations. 

In the forthcoming derivations the estimate (6.5) will be used to bound the probability of a 
single bad event but, more importantly, to decouple various bad events. However, it will not be 
always advantageous to estimate ZL^f3{Ak) directly — often we will have to further decompose Ak 
into smaller events. Then we will use the well-known subadditivity property: 

Lemma 6.3 (Subadditivity) Consider the events A and {Ak)keJ^ that depend only on configu- 
rations in a box F fl where V is as specified above. If Ac Uite^ "^k' ^^^^ 



ZlM) < ( E • (6.6) 



Proof The claim is (presumably) standard; we provide a short proof for reader's convenience. 
Clearly, Zl^p{A) is equal to the Z^^^^-multiple of the expectation of HteT '^^tiA)- Now, using 
the bound < J2keJ^'^MAk) get 

where the collection {kt)^^j G provides the assignment of a /ct G ^ for each of the 
translates of the .4-events and where the sum is over all such assignments. Applying Theo- 
rem 6.2, Zij3 times the expectation on the right-hand side of (6.7) is bounded by the product of 

Zi.piAkt)^^^'^^ over all t G T. But then each kt can be independently summed over whereby the 
desired relation (6.6) follows. □ 

6.2 Distinct types of badness. 

The estimate of the probability of bad events — defined right after Definition 1 in Sect. 3.1 — will 
require partitioning this event into further categories. A priori, we will distinguish two types of 
badness according to which violation of the aforementioned conditions in Definition 1 is high- 
hghted. Specifically, we define the events 

Be = {S: \sj:"'^ - S^'^^J > r for some r,r + e„ € Ab} (6.8) 

and 

Bsw = B\ Be. (6.9) 
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Here is the cube of {B + 1)^ sites with the "lowest left-most" site at the origin (this is where 
the prototype bad event B was defined). The idea behind this spUtting is that for the configurations 
in Be there is an energetic "disaster" while for those in Bsw the spin-wave approximation is still 
good but we are not particularly close to any free-energy minimum. 

Unfortunately, the event Bsw is still too complex to be estimated directly because after dissem- 
ination all over the torus, the resulting partition function ^l,/3('Ssw) does not end up being the 
type featured in Sects. 4-5. This is directly related to the fact — eluded to earlier — that there are a 
myriad of ground states in these models. Thus we will have to work a bit in order to parcel Bsw 
into events which after eventual dissemination over the torus lead to partition functions of a type 
discussed in the previous sections. 

In order to motivate the forthcoming definitions, let us categorize, somewhat more precisely 
than in Sect. 1.3, the ground states of the model (1.1). To avoid intricacies due to boundary 
conditions, we will restrict ourselves to toroidal geometry. First, all constant spin configurations 
minimize the energy. Second, more ground states can be generated from a homogeneous config- 
uration by picking a lattice direction, a, and a sequence of planes orthogonal to e^, and reflecting 
all spins in these planes through a-th of the vectors a,b,c. These statements are more or less 
fully justified by Proposition 6.5 below. Of course, when we split ;Ssw into the homogeneous and 
inhomogeneous parts, we will not try to keep track of all planes of reflection that can occur — one 

(i) 

will be sufficient (this is the basis of the event B\,'- below). 

The decomposition of the event Bsw will involve all of our basic scales: For a given k > 0, 
r > and 5 > we let 

A = 1?^. (6.10) 

K 

(As we will see in Sect. 6.3, this will be the A for which we will use the results of spin-wave 

analysis from Sect. 4.) Let us parametrize the spins using the angle variables 9r. Fix an integer 
s > 1 and \et9^, . . . ,6g be s points uniformly spaced on the unit circle. The first part we want to 

identify from ;Ssw are the nearly homogenous configurations: For each i G {1, . . . , s}, let Bq^ 
denote the event that the block A^ is bad and that \0r — 0*\ < A holds for all r G Ag. 

The complementary part of Bsw will feature a particular kind of inhomogeneity: Fix an 
a G {1, 2, 3} and let j G {1,2,..., B} and let Mj denote the plane in where all r have 
the a-th component equal to j. Fix i G {1, . . . , s}. If the angle between 0i and the a-th of the 
vectors a, b, c is within (— k, /t) or (vr — k, vr + k), then we set B^^j = 0. For the other i we let B^^^j 
denote the set of all configurations S G Bsw such that \9r — 0*\ < A holds for all r G A^ n EIj_i 
and \9r - 6*\ < A for all r G Ag n M^. Here 6* denotes the angle 0* reflected through the a-th 
of the vectors a, b and c. 

Remark 7. Let us reiterate that, by definition, we have Bq^ C B for all i and B^^^j = for all i 
whose 6i is too "near" the a-th of the vectors a, b, c. These facts will be useful when we estimate 
the associated partition functions Z^jj^Bq^) and Zig{B^^-) in Sect. 6.3. 

It remains to show that the union of these events contains all of ^Ssw^ 



LONG-RANGE ORDER IN 120° -MODEL 



29 



Theorem 6.4 Consider the 120° -model and let the events BsW' and Bj- be as defined 
above. Suppose that V, k, B and s are such that By/V <^ k <^ 1 and sA > Air. Then 

BswCU(^?^U U U^S)- (6.11) 

i=l a=l,2,3 j=l 

Remark 8. In the above and in what is to follow (and to a certain extent retroactively) we employ 
the symbol "<C" in our hypotheses according to the standard fashion: "if a <C b. .." means "if 
the ratio a/6 is bounded by a sufficiently small numerical constant which is uniform in any of the 
other parameters mentioned " 

Remark 9. The inclusion (6.1 1) justifies our previous claim that the only spin-wave calculations 
we need to do are against either homogeneous or stratified background. Indeed, by Lemma 6.3, 
to estimate the probability of Bsw we will only need to estimate the constrained partition func- 
tions Zl^p{Bq^) and ZL^p{B'^^j). The former leads directly to homogeneous spin- wave calcu- 
lations from Sect. 4.1; the latter will require further disemination of the pair of {j — 1, j)-th 
planes in the a-direction which results in exactly the stratified background configuration treated 
in Sect. 4.2. See Lemmas 6.9 and 6.10 for details. 

The proof of Theorem 6.4 commences by considering an elementary cube in T^, say K = 
{0, 1}'^, and classifying all spin configurations on K that are "nearly" a ground state but which 
are not near any of the six "priviledged" directions w^; see (2.1) for the corresponding definition. 
The precise statement is as follows: 

Proposition 6.5 Let T and k be such that ^/T <C k ^ 1. Let 6 = (Or) be a configuration of 
angle variables on K such that the corresponding spins Sr satisfy the energy constraints (3.1) /or 
all pairs of nearest neighbors on IK but such that not all of the spins are within angle k of one 
particular w,-. Letr G K. Then (exactly) one of the following is true: 

(1) {Or' - Orl < AF/k for all r' G K 

(2) There exists an a G {1,2,3} such that \6r' — 0r\ < 4:T/k holds for all r' G K with 
r — r' 1. Ga, while for the remaining r G K we have {Or' — 9r\ < / k, where Or is obtained 
from Or by reflection through the a-th of the vectors a, b, c. 

Remark 10. Setting F = (and k = 0) in this statement justifies Fig. 1, which shows four exam- 
ples of ground state configurations on an elementary cube. The reason why we explicitly exclude 
the "almost" constant configurations which point near one of w^ is that, for these situations, the 
energy constraint would permit fluctuations that are of order v^. 

The proof of Proposition 6.5 will involve a couple of lemmas. First, let us characterize the 
consequence of the energy constraint (3.1) for a single bond: 

Lemma 6.6 Let a G {1, 2, 3} and consider a nearest-neighbor bond {r,r') parallel to Sq. Let Or 
and Or' be two angle variables such that the corresponding spins satisfy IS^"^ — S^f ^ | < F. Then 
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either \dr — 6ri\ < 7r-^r/2 or \9r — 0r'\ < 7ryT72, where Or' is obtained from Or' by reflection 
through the a-th of the vectors a, b, c. 

Proof. Without loss of generality, we will assume that a = I. Now, if 6, 9' G [0, vr] are two 
angles with \6 — 0'\ = e, then the trig identity | cos 6 — cos 0'\ = 2| sin(^^)| sin(^i^) and some 
optimization show that 

\cos9-cose'\ > 2sin2(§) > 2eV7^^ (6.12) 

where we used that e/2 G [0,7r/2]. But the left hand side is exactly |S'/^'* — S^j^\ which by 
assumption is less than T. A simple algebra now shows that then e. = \9 — 6'\ < ti^Y 12. This 
proves the claim in the case when both Q and Q' have the same sign; the opposite case is handled 
by reflection through the x axis. □ 

Next we will extend this to a similar control of lattice plaquettes: 

Lemma 6.7 Let 'Lbe a lattice plaquette and let C = 67r/ \/2. Let F <C 1 and let Or' , r' G L 
denote angle variables such that the energy constraint (3.1) holds for all four bonds. Then, for 

any particular r G L, either all 6^', r' G L, are within C\IV from Or or one neighbor ofr satisfies 
this condition while, on the other side of the plaquette, the other two spins are within C\/V from 
the corresponding reflection of Or. 

Proof. The proof is based on Lemma 6.6. To make the reference to this lemma easier, let us say 
that a reflection occurs for the pair (r , r') if the latter possibility in Lemma 6.6 applies. Let L be a 
lattice plaquette. Since permutations of coordinate directions can be matched with permuting the 

roles of a,b, c, we can as well assume that L is an .xy-plaquette, i.e., L = {r.r + ex-,r + ty^r + 
Cx- + Cy}. The analysis proceeds by checking various cases of increasing level of complexity. To 
simplify the formulas, let us abbreviate the error constant from Lemma 6.6 by C = 7ryT/2. 

CASE 1 : No reflection occurs for both of the bonds emanating from r. Lemma 6.6 then implies 
that both ^r+Sa: ^r+g;, within C, from Or- Now if a reflection does not occur on either of 

the two remaining bonds, then the spin at r + e^- + e^ is within 2^ of Or and we are done. The 
remaining possibility would be a reflection on both of these bonds. But then ^r+ea;+ej, is within 2(^ 
of both —Or and {—\ — Or) which is impossible once AC, < 

CASE 2: Reflection occurs for exactly one bond emanating from r, say the horizontal bond 
from r. The only case we need to consider is when reflection occurs for the "other" vertical 
bond and does not for the "other" horizontal bond. But then Or+i^+iy is within 2^ of both Or 

and {—'^ + 0r) which is again impossible once 4^ < 

CASE 3: Reflection occurs for both bonds emanating from r. Clearly, following the path through 
r+Qx tell us that Or+e^+^y is within 2^ of either —Or or {—^+0r) while the passage through r+Qy 
tells us that Or+s^-\-e:y is within 2^ of ±(^ + 0^). Checking the cases shows that, \£ AC, < ^, 
this is only possible when reflection occurs for all bonds around the plaquette and when Or is 
within 2^ of one of the angles 0, ^, vr or Let us check the case when Or ~ 0. Then Or-\.^^ is 
within 5C of Or and both ^r+S,, ^nd ^r+Sx+S,; within 6^ from (— ^ — Or). A similar argument 
handles the remaining cases. 
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Inspecting the above derivations, we see that the worst-case fluctuation from one of the two 
situations described in the statement of the lemma is by 6^ = C\/r. This finishes the proof. □ 

Now we are ready to characterize the "near" ground states on elementary cubes: 

Proof of Proposition 6.5. Lemma 6.7 immediately implies that any configuration S satisfying 
the energy constraints (3.1) on K is one of the types featured in the statement of the proposition 
(resp., Fig. 1) to within errors C"\/r for some numerical constant C' . Indeed, either all angle 
variables are within C"\/f of some particular angle or not. If not, then there must be a pair of 
nearest neighbors {r,r'), say parallel to ei, where a reflection occurred. Then 9r' is within CVt 
of —dr. Moreover, choosing C' ^ C allows us to assume that both j^^l and Ivr— exceed 2C\/r 
and thus both plaquettes in K containing the bond {r,r') will have to be of a "mixed" type. But, 
again by Lemma 6.7, the two perpendicular yz-plaquettes cannot be of a "mixed" type. This 
imphes the characterization in the statement of the proposition with the errors bounded by C'Vr. 

It remains to show that the errors are in fact only proportional to T (cf Remark 10). Here 
we will use the following refinement of Lemma 6.6: If 9, 9' € [0, vr) satisfy the energy con- 
straint I cos 9 — cos 9'\ <T but are not within angle k ^ 1 of the x-ground states, then 

\9-9'\<r/K. (6.13) 

Indeed, the Mean Value Theorem gives us that | cos 9 — cos 9'\ = 2 sm{9"){9 — 9') where 9" hes 
between 9 and 9'. Hence sin(0") > sinK which for k <C 1 exceeds k/2. Using that \ cos9 — 
cos 9'\ < r, the bound (6.13) directly follows. 

The improved error bound is now a simple consequence of (6.13). Let us first consider the 
"nearly" homogeneous situations. Since all angle variables are to be away from the ground state, 
(6.13) implies that for each bond the 6'r's will differ only by at most T / n. Hence, all 9r& on the 
cube must be within SF/k of one of them which proves the claim in this case. The "mixed" con- 
figurations whose both types point away from any of the ground states are handled analogously, 
so we only have to consider the case when each type is within angle 2 k of a different ground 
state. A generic situation of this kind is when the "bottom" xy-plaquette of K is occupied by a 
configuration 9 k, {) while the "top" xy-plaquette is occupied by a configuration 9 k, Then 
the observation (6.13) constrains the size of the fluctuations to less than V/k along the following 
bonds: The e^^-bonds in the "bottom" xy-plaquette, all of the vertical bonds and the e^-bonds in 
the "top" xy-plaquette. It is easy to check that, from r, one can get to all sites of IC in at most 
four steps, so all 9ri are within less than 4r/K of 9r or the corresponding reflection. □ 

Finally, we are ready to prove Theorem 6.4: 

Proof of Theorem 6.4. Consider a spin configuration S on such that i3sw occurs and let 9r be 
the corresponding angle variables. Suppose first that one of the 0r's makes an angle at least 2k 
with all of the w,-, r = 1, . . . , 6. Applying Lemma 6.6 along with the fact that B\/T ^ k, we 
find out that all 9/s will be make an angle at least k with any of the w,-. Proposition 6.5 then 
guarantees that any elementary cube has a layered structure with the 9/8 more or less constant 
in both layers. Since the maximal fluctuation in each elementary cube is at most 4r/K, it is not 
more than SB-times that — ^i.e., A in (6.10) — for any pair of spins in Ag. 
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Now the bound sA > An ensures that the consecutive 9f (which we used to define the events 
Bq^ and B^^j) are within less than A/2 from each other. Thus, if all spins point in about the same 

direction they must all be within A of some Of — which implies that S G Bq^ — or there are two 
consecutive layers, say j — I and j, in the a-th lattice direction where a reflection from 6* to 6* 

(i) 

occurs. In the latter case we have S £ ^a j- "^^^^ proves (6.11) for those S G Bsw for which 
at least one of the spins is farther than 2k (in the angular distance) from any of the six preferred 
directions w,-. 

It remains to deal with the situations in which all spins are within 2k of some w,- (possibly 
different r for different spins). Clearly, the latter cannot be the same for all spins because of 
the inclusion Bsw C B, and so there must be a pair of spins where the type of ground state is 
different at the endpoints. But then we can still use Proposition 6.5 for the elementary cubes 
containing this bond, and then the cubes next to these and so on. In this way we conclude that 
the endpoints of this bond belong to two parallel planes of sites where the spins do not fluctuate 
by more than 2B times AT/k about a single direction in one plane and its reflection in the other. 

Hence S belongs to one ofthe^J^s. □ 
6.3 Proof of Theorem 3.1. 

We begin with an estimate of the partition function for event Be- 

Lemma 6.8 Let k > be fixed. There exist constants C4 G (0, 00) and (5 > such that if PJ, 
A = \2BV /Kand 5 satisfy the bounds (4.2) then 

limsup < B\c,(5Jr/^e--^P^^\ (6.14) 

Proof. We will derive an upper bound Z^^j^iB^) and a lower bound on Z^^p. The former is 
essentially an immediate consequence of the definition of B^. Indeed, on B^e at least one of the 
pairs of nearest neighbors in contributes at least i(/?J)r^ to the total energy. Thus, after 
dissemination of B^. all over the torus, the spin configurations are constrained to satisfy 

/?^L(Sf)>^(/?J)r2(^)'. (615) 

It follows that 

^l,/3(He)(^/^)' < QB\2TTf'e-'^^J^\ (6.16) 

where the factor QB^ bounds the number of places where the "excited" bond can occur within 

and (27r)^^ is the total "phase volume" of all configurations in A^. 

Next we need to derive a lower bound on Z^^p. Here we will write the partition function as 
an integral of e~^-^^ ; a lower bound can then be obtained by inserting the indicator that all angle 
variables are within A of 0°. This yields 

> y— j e (6.17) 



LONG-RANGE ORDER IN 120° -MODEL 33 

where F^^'*(0°) is the quantity from (4.1). Choosing e > and letting 5 be such that Theorem 4. 1 
holds, we thus get 

liminf(Z,,,)V^^>(|^)'^%-^(°°)-, (6.18) 

where F denotes the spin-wave free energy (2.11). Combining (6.15) and (6.18) and letting C4 
absorb all factors independent of B and (3 J, the desired bound (6.14) is proved. □ 

Remark 11. Since the event depends only on V, the appearance of k in the assumptions of 
Lemma 6.8 may seem unnecessary. However, some conditions on F, B and (3 J are still needed 
to derive the lower bound in (6.18) and the advantage of the present form is that now all lemmas 
in this section are proved under more or less the same assumptions. 

Next we will attend to the event 6sw- In light of Theorem 6.4 and Lemma 6.3, we can focus 
directly on the events B^^ and B^^-. We will begin with the former of the two: 

Lemma 6.9 Let k > be fixed. There exist numbers p\ (k) > and > such that if (3 J 
and A = 12BT / n satisfy the bounds (4.2) with this 5 and if BY <^ k <^1, then 



limsup ( ^h^llV"-" < e-Pi(«)^ z = 1, . . . , (6.19) 



Proof. To summarize the situation, on Bq \ all angle variables 9r in the block Ab are within A 

of 9*. If we now consider the multiply reflected event associated with Bq \ the same will be true 
about all spins on T^,. Let e > and let 5 > be as in Theorem 4.1. Then 

limsupZ,,,(i3«)V^^ < (^f'e-m)^^. (6.20) 

Using (6.18) we thus conclude 



L-+0O ^ ZL,f3 

(i) 

It remains to adjust e so that the exponent is negative. Here we first note that Bq is empty un- 
less 9f is at least k away from any of the ground state (indeed, otherwise the configuration fails 
to be in B, which by definition contains Bq"^). Applying Corollary 5.2, F{9*) exceeds F(0°) by a 
uniformly positive amount, denoted by 2pi(k), whenever 9* is at least k away from the minimiz- 
ing angles. Now choose e < ^Pi{k) and let 6 be the corresponding quantity from Theorem 4.1. 
Then the right-hand side of (6.21) is indeed less than e~^^^'^\ proving the desired claim. □ 

(i) 

Similarly, we have to derive a corresponding bound for the events B)^ ■: 
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Lemma 6.10 Let k > be fixed. There exist numbers P2{n) > and S > such that if (3 J 
and A = 12BT/k satisfy the bound (4.26) with this S and ifBT <C k <C 1, then 

limsup(^^#S)y/^^<,-p.(.)/B (6.22) 
holds for all a G {1, 2, 3}, all j € {1, 2, . . . , B} and all i G {1, . . . , s}. 

Remark 12. We assure the reader that the 1/B in the exponent is no cause for alarm; in accord 
with (6.6), the relevant object from Lemma 6.10 is the right-hand side raised to power B^. 

(i) 

Proof of Lemma 6.10. Recall that, on B Jy all Or for r in the plane H ]HIj_i are within 
a constant times BT/k of 0^, while those in the neighboring plane n Hj are within the 
same distance of the reflected angle 0*. After dissemination over the torus, which is what gives 
rise to the quantity Zl^p{b'^^j), the same will be true about the spins in the entire planes Mj^i, 
resp., H^ , and also about their translates by integer multiples of B in the orthogonal direction. 
However, we cannot yet use the spin-wave calculation; instead, we have to use Theorem 6.2 again 
to disseminate the two-plane alternating pattern all over the torus. This yields 



where b''^),- is the event in A^ that the Qr are within A of ^* in even translates of IHIj_i and of Q\ 
in odd translates of BIj_i. 

Now the partition function can be estimated using Theorem 4.5 and we thus get 

limsup i^b^^^i^ < e-[^(er)-F(o=)-2.]^ (6_24) 

But Theorem 5.5 shows that Fa{0*) — F{0°) > C3 > for some C3 = c^{k) for all i for which 9* 
is at least K-away from any of the minimizing angles associated with "stratification" direction a, 
while, by definition, B^;^\ = for those i that fail this condition. Hence, if we choose e > so 
small that p2{n) = 2(c3 — 2e) > and let S be the corresponding constant from Theorem 4.5, 
then (6.23-6.24) imply (6.22) as desired. □ 

Proof of Theorem 3.1. Let k > and let ^ > be the minimum of the corresponding num- 
bers from Lemmas 6.8-6.10. Fix an r/ G (0, 1). We claim that for each sufficiently large (5, 
there exist numbers B and T such that the bounds (4.2) and (4.26) for A = 12BT / n hold, the 
inequality By/V <C n can be achieved and the bound 

B\c,pj)''''\-lP'^' + |:e-^^''i('^) + ^e-^^^^('') < (6.25) 

is true. Indeed, we can for instance take B = log /3 and F = /3~ 12 and note that, for these 
choices, the left-hand side will eventually decrease with /3. 
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Now choose s such that sA > An but sA < Stt. Then the definitions (6.8-6.9) of events Be 
and Hsw> the decomposition of ;Bsw from Theorem 6.4, the chessboard estimate and the (subad- 
ditivity) Lemma 6.6 imply that ^L,/3i'&ti (B) PI • • • n -dtm i^)) will be bounded by rj'^, where 



V Zr. R / V Zt, R / 



E E 



=1 a=l,2,3 i=l 



1 . (6.26) 



By Lemmas 6.8-6.10, the fact that s < Sir/ A and (6.25), it follows that limsup^^^o rn < f]. 
Hence there exists a number Lq G (0, oo) such that riL < rj for all L > Lq. But for L > Lq, 
the probability PL,/3(i9ti {B) PI ■ ■ ■ fl {B)) is bounded by rf^ uniformly in m and the choice of 
the vectors ti , . . . , tm- This proves the desired claim and thus also finishes the proof of our main 
result (Theorem 2.1). □ 



7. Spherical models 

Here we present the proof that the spherical version of the 120° -model has no phase transition 
at any positive temperature. This demonstrates the failure of the naive spin-wave arguments and, 
particularly, the infrared bounds. 

Spherical models, very popular in the 1950-60, were conceived of by Berlin-Kac [2] as conve- 
nient approximations of the statistical mechanical systems which are more amenable to explicit 
computations. (On the mathematics side, the topic received a new wave of interest in the 1980's 
through the rigorous versions of 1/n expansion.) To construct a spherical version of a given spin 
system, we use the same Hamiltonian but ascribe different meaning to the spin variables. In par- 
ticular, the local a priori constraints on the spin variables are relaxed and are replaced by a global 
constraint. For instance, for the Ising model with Hamiltonian = — ^, (7^(7^' we have 
(Tr = ±1 and thus cr^ = 1 for all r. The spherical version has the same interaction Hamiltonian 
but now we only require that (l/N) = 1, where N denotes the total number of spins. 

Often enough, these models are further simplified by stipulating that the constraint only needs 
to be satisfied in the mean and may thus be enforced by Langrange multipliers. The latter type is 
often referred to as the mean spherical model. This version usually turns out to be pretty much 
the same in most aspects, see [27] for some discussion. Here we will go the mean-spherical route 
partially because the resulting analysis is simpler, but also because the analogy to pure spin-wave 
theory is more pronounced in this case. We refer to [44, Section 11.11] for more references and 
further discussion. 

Thus, we will take (1.10) as our basic Hamiltonian along with an additional term to enforce 
the required constraints. However now it is understood that the spin variables are no-longer con- 
strained to the unit circle; the integration takes place over all of M^. The constraining term reads 
—jj, J2r a ("S'r"^ ) ^ but uow (unfortunatcly) 5^"^ refers to the Cartesian component of the spin. This 
means that we will have to rewrite the Hamiltonian in terms of the Cartesian components of S. 
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The key to the mean-spherical approximation is that for arbitrary > the partition function 
can be solved exactly by translations to spin- wave variables. Then jj, is supposed to be adjusted 
so that the relevant constraint is enforced. As we shall see, if there is an infrared divergence, 
this adjustment is easy and everything is analytic in (3. In the opposite case, there may be a 
condensation at large /3 and if so, one may conclude — with a lot of apologies — that a phase 
transition has occurred. The primary conclusion of this section is that the latter possibility does 
not materiaUze in the model at hand. 

Now we are ready to describe the spherical version of the 120° -model. The Hamiltonian on 
torus is given by 



+ 



(7.1) 



(x) (v) 

where Sr and Sr are now unrestricted real variables a priori distributed according to the 



Lebesgue measure on R. The constraint is represented by the quantity 



(7.2) 



The associated Gibbs measure is given in terms of the Radon-Nikodym derivative with respect to 
the Lebesgue measure on (R^)^^, which is simply a properly normalized erl^-^^~>^-^^ . We will 
denote the expectation with respect to the resulting thermal state by {— ) l,^,/*- 

Theorem 7.1 Consider the spherical 120P-model with the Hamiltonian (7.1) and let {—)L,/3,jj, 
denote the corresponding thermal state for the chemical potential p. Then there exists a positive, 
real-analytic function /x* : [0, oo) — > (0, oo) such that for each /3 G (0, oo) — and ji set to //*(/?) — 
the following is true: The constraint is satisfied on average, 



1 



there is no long range order. 



lim 

L—*oo 



= 0, 



(7.3) 



(7.4) 



and the limiting measure exhibits a clustering property. 



(7.5) 



for any a, a' G {x, y}. Moreover, the limiting free energy is a (real) analytic function of f3. 

Proof. As usual, our first goal will be to calculate the limiting free energy as a function of (3 
and /Lt. Let Zl{P,p) denote the integral of g-|3•^L-^^■^L ^vith respect to the Lebesgue measure 
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on (M^)^^. In order to compute Zl{P, //) we transform to the Fourier modes in which case the 
spin- wave Hamiltonian (including the constraint) is seen to be given by 

+ \si'^\' [f (^2 + £^3) + A] + ^{si^^S^^'l + S'flsi^^) [E2 -Es]y (7.6) 

Here S^^/^ is just the complex conjugate of sj^"\ the symbol Ea abbreviates the usual Ea{k) = 

|1 — e'^^^P and A is definedby /3JA/2 = ^. In terms of the two-component variable (4t 

the right-hand side (without the /3J/2 prefactor, of course) can be written as a quadratic form 

with matrix A1 -|- ©(A:), where 

q,j^.(e, + \{E2 + E^) f(E2-^3A .77. 

^{E2-E^) \{E2 + E^) )' 

In this notation the integrals are readily performed with the limiting free energy F{I3, A) — which 
to within a sign is the hmit lim^^oo log [3JX/2) — given by 

no, A) = log 1^ 4 1^ log det [A1 + e(t)] . (7.8) 

Here the integral converges as long as A > 0. 

Our next goal is to find the function ji^, for which (7.3) holds. Using standard relation between 
free energy and expectation, the constraint equation becomes 

where 

5sp(A:) = lim {\Sk?)L,p,^. = (/3J)-^Tr[A1 + Q{k)]-^ (7.10) 

is the so-called structure factor. As long as A > 0, the derivative -^F(P, A) is finite and inde- 
pendent of P and thus (7.9) defines a function A i-^ /3^(A). A moment's thought shows that this 
function is strictly decreasing and hence locally invertible. However, before we plug the inverse 
back into (7.8), we need to estabUsh the range of values that /?*(A) can take. In particular, we ask 
whether /3*(A) diverges as A J, 0. 

Examining the constraint equation in detail, the crucial issue boils down to convergence/diver- 
gence of the momentum-space integral of the structure factor 

osp(k) oc . (7.11) 

E1E2 + EiEs + E2ES 

It turns out that the integral of SsY>{k) diverges although this is not apparent by naive power 
counting. Indeed, the primary source of the divergence is not the origin but the coordinate axes. 
This is seen by an easy lower bound on Ssp{k): Fix ks to a non-zero number and note that we 
can discard the Ei and E2 from the numerator. Second, the term E1E2 in the denominator is 
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while 



bounded above by a constant times Ei + E2. Hence, the calculations boil down to the integral of 
(£'1 + E2)~^ with respect to k2 and kz, which is manifestly divergent. 

The above reasoning shows that A ^ /3*(A) takes all positive real values as A sweeps through 
the positive real line and hence the inverse f3 i-^- A*(/3) is defined for all (5 G [0, 00). Moreover, 
for A > the function A 1— > ( A) is analytic in a small neighborhood of the real fine and hence so 
\s (3 ^ X~k{(3). The desired function then arises by setting = l3J\i,{(3)/2, which satisfies 

(7.3) by construction. Furthermore, plugging A^(/3) for A in F{\, (3) proves that the free energy 
is real analytic in /3. In order to prove also (7.4-7.5), we just need to note that (7.6) implies that 

the correlator (^i^^^^J"'^) l., is exactly the (a, a')-th matrix element of (/3 J) ^ [A1 + 0(^)] ^. 
But then 

which by the Riemann-Lebesgue lemma and the fact that X-k{l3) is strictly positive for any P G 
[0, cxd) tends to zero as |r — r'| — 00. □ 

Remark 13. The last expression of the proof indicates that the correlations decay (at least) expo- 
nentially fast. However, as is seen from (7.11), the angular dependence of the resulting correlation 
length is fairly complicated. In particular, there may be directions in which the quadratic approx- 
imation of &{k) vanishes in which case more than one pole in the "complex \k\" plane (instead 
of the usual single pole) jointly contribute to the integral. 

We conclude with a remark concerning the relation of these findings to the actual systems of 
interest. For the spherical model, the so called structure factor Ssp{k) = {\Sk\'^)i3^^ can explicitly 
be computed, cf (7.10). As was established in [22,23,25] for a general class of nearest-neighbor 
ferromagnetic systems (including the one discussed in the present work) the spherical rendition 
of the structure factor with ^ = provides a bound on the structure factor SA{k) (namely, the 
two-point correlation function in A: -representation) of the actual system, 

5a(*)<<Ssp(*)|^=o- C^-l^) 

This is the basis of the infrared-bound technology which uses the convergence of the integrated 
bound to estabUsh long-range order at low temperatures. 

Here, the low momentum behavior of the spherical structure factor together with the rigorous 
as well as non-rigorous results relating Ssp to 5a (including in particular [28]) strongly suggest 
a disordering due to long wave-length fluctuations. It is usually the case that these are reUable 
indicators for the behavior of the actual system. Evidently, as the results of this work show, the 
present cases are exceptional. 
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